Polarized superfiuids near their tricritical point 



o 
o 

(N 

X) 

(N 



O 

o 
I 



O 

> 

(N 

0^ 
p 

O 
00 

O 



X 




Daniel E. Sheehy 
Department of Physics and Astronomy, Louisiana State University, Baton Rouge, Louisiana 70803 

(Dated: July 6, 2008; Revised Dec. 30, 2008) 

We obtain a variety of predictions for the properties of population-imbalanced (or polarized) 
fermionic superfiuids near their tricritical point. In the vicinity of the high-symmetry tricritical 
point, observable quantities such as the cloud shape, heat capacity, local polarization and correlation 
length should exhibit distinct behavior arising from the tricritical scaling laws, as well as logarithmic 
corrections to scaling reflecting the marginal nature of interactions. 



I. INTRODUCTION 

One of the most exciting recent developments in atomic 
physics has been the achievement ofpaired superfluidity 
of fermionic atomic gases [l|, 0, B 0, a B 0] ■ Such super- 
fluidity, arising from attractive interactions between two 
fcrmion species mediated by a magnetic-field tuned Fesh- 
bach resonance, can be continously experimentally tuned 
(by adjusting the external magnetic field) from the BCS 
limit of weak pairing to the BEC limit of strong pair- 
ingii. 

In fact, the smooth crossover between these two lim- 
its only occurs for an equal number of the two fermion 
species and any population imbalance (or polarization) 
interrupts the BEC-BCS crossover, as seen in recent ex- 
periments [S [nl [13, [H- Theoretical work [13, [H, 
[iq . \VA [isl [ig . [2q | on such polarized Fermi gases pre- 
dicts a rich ground-state phase diagram including an 
FFLO phase [22| in the positive detuning BCS regime, 
a strongly-interacting polarized normal Fermi gas at 
large polarization [23|, a regime of magnetic superflu- 
idity (consisting of tightly-bound molecules coexisting 
with a single-spin Fermi sea, analogous to '^He-'^He mix- 
tures ^, 25]), and a regime of phase separation. 

Although a consistent picture of the ground-state 
phase diagram has arisen from mean-field theory [l^, [l8| 
and Monte Carlo i20i] r esults, the finite-temperature 
phase diagram [i3j[ll,[3 is less well understood, partic- 
ularly near the unitary point where the interspecies scat- 
tering length fls diverges. For any detuning, at sufRcently 
low polarization or chemical potential difference H we ex- 
pect a finite-temperature phase transition at which the 
superfluid order parameter vanishes continuously. Con- 
versely, at low temperature, superfluidity is destroyed in 
a flrst-order fashion [26[ with increasing H. Across this 
phase transition at fixed H, the polarization or magne- 
tization M = nij — ni (with Ua- the density of spin-cr 
fermions) jumps discontinously; at fixed M one finds a 
regime of phase separation [l^, [1^1 , as seen in Fig. [TId. 
It is natural to suppose that these first and second-order 
regimes are connected by a tricritical point [2l| (TP), as 
which occurs in BCS superconductors under an imposed 
magnetic field and in ^He-'*He mixtures. 

Indeed, recent experiments [l^ ^^ have reported ev- 
idence for a TP in the unitary regime, so that an ex- 
ploration of additional tricritical phenomenology is of 
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FIG. 1: (Color Online) (a) Phase diagram vs. chemical po- 
tential difference H and temperature T near the tricritical 
point (TP) at Ht,Tt (solid point), showing magnetic super- 
fluid (SFm) and polarized Fermi liquid (N) regions separated 
by a continuous transition at T > Tt (thin line) and a first- 
order transition at T < Tt (thick line), (b) Phase diagram 
vs. magnetization M and T, near the TP at Mt,Tt, show- 
ing the coexistence region (shaded; bounded by Md and 
Mc2), with the coexistence curve asymptotically satisfying 
M,2- 
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A/tln|i| with t={T- Tt)/Tt. 



considerable interest. Narrow-resonance models 18[ pre- 
dict that a line of tricritical points crosses the phase 
diagram as a function of Feshbach resonance detuning, 
terminatin g at a quantum tricritical point in the deep 
BEC limit [T^. However, in the experimentally-relevant 
wide-resonance limit, it is possible that strong correla- 
tions may interrupt the tricritical point, as found in the 
Blume- Emery-Griffiths model [2J]. Thus, our aim is to 
devise novel experimenal signatures of a tricritical point, 
to aid in establishing the phase diagram of polarized 
Fermi gases. 

We thus proceed by assuming that the phase diagram 
of polarized superfiuids possesses a TP. If so, then, quite 
generally, the phase diagram near such a tricritical point 
will resemble Fig.[Tl neglecting the possibility of an FFLO 
state, which, at least within mean-field theory, is re- 
stricted to a thin window of ff or M values [la, llMl 
in the BCS regime. Having made this assumption, we 
shall make predictions for the behavior, near the TP, of 
various observable quantities in cold atom experiments, 
such as the local polarization (magnetization) , molecular 
or pair density, and the heat capacity. Our results are 
based on the analysis of a sixth-order Ginzburg-Landau 



(GL) free energy via mean-field theory and the renormal- 
ization group (RG) and are expected to apply generally 
to finite- temperature TPs in polarized Fermi gases [27| . 
This manuscript is organized as follows. In Sec. |TT1 we 
introduce the GL model of tricritical points and discuss 
our principal results. In Sec. lIIIi we derive the coefficients 
of the GL model within a mean-field analysis of the one- 
channel model of resonantly-interacting polarized Fermi 
gases. In Sec. |lVl we derive the RG equations for the 
sixth-order GL model. In Sec. |V] we use the RG to de- 
rive equations for the variation of the superfluid order 
parameter below the phase transition in the vicinity of 
the tricritical point and show how these results imply a 
characteristic steep cloud shape in a trapped polarized 
gas. In Sec. I VII we derive a prediction for the jump in 
the magnetization M across the first-order phase bound- 
ary near the tricritical point, a quantity that translates, 
in a trapped polarized superfluid, to a jump in M as a 
function of radius. In Sec. IVIII we derive equations for 
the heat capacity above and below the transition near 
the tricritical point. In Sec. IVIlIl we describe the diver- 
gence of the order parameter correlation length near the 
tricritical point before concluding in Sec. IIXI 



II. MODEL AND PRINCIPAL RESULTS 

As we have noted, a tricritical point in the phase di- 
agram of polarized Fermi gases separates first-order and 
continuous phase transitions of the superfluid order pa- 
rameter ■)/), with \'i}j\^ essentially representing the density 
of condensed molecular pairs [2q. The behavior near 
the tricritical point can be captured using the following 
sixth-order free-energy functional: 
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+ 2^-1^1' + 4"i^r+g«iV'r , (1) 



where mi, = 2m is the molecular mass, and r, u and v are 
T (temperature) and H (chemical potential difference) 
dependent coefficients. 

Below, we show how Eq. ^ can be obtained within 
a mean-field analysis of the standard one-channel model 
of resonantly-interacting Fermi gases, allowing us to de- 
rive mean-field predictions for the coefficients r, u, v. 
Although mean-field theory breaks down in the unitary 
regime, we nonetheless expect Eq. ([T]) to correctly de- 
scribe the vicinity of the tricritical point, but with un- 
known coefficients. In the present section, we proceed by 
reviewing the mean-field phase diagram of Eq. ([T]). 

The basic tricritical phenomenology (including the 
phase diagram in the u — r plane, Fig. [2]) follows from 
analyzing Eq. ([1]) in the mean-field approximation, by 
minimizing with respect to an assumed spatially-uniform 
-0. With a spatially uniform -0, the mean-field free energy 
is: 



F=\rW + \uW^\v\^\^ 
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FIG. 2: (Color Online) Phase diagram of F (with v = 
1), Eq. lU, showing superfluid and normal (nonsuperfluid) 
phases separated by a first-order boundary for u < and a 
continuous boundary for it > with the tricritical point at 
the origin. 

where we note that, in going from Eq. ^ to Eq. ([2]), we 
have clearly taken the system volume to be unity. Along 
with the trivial stationarity condition = 0, one finds 
the nontrivial stationarity condition 



O = r + u|0p + u|0|^ 



that can be solved to yield 



101^ = -— \\/v? ~ 4rv - 
2v 



(3) 



(4) 



(2) 



We begin with the regime, u > 0, in which the phase 
transition is continuous, that occurs for low H . Since we 
expect V > (necessary to stabilize the tricritical point), 
we see that Eq. (|4]) only represents a physical solution 
for r < 0, so that the stable solution is V" = for r > 0. 
This nonsuperfiuid, or normal, regime is characterized 
by a vanishing pair (or molecular) density n™ = |0p. 
For r < 0, in the u > regime superfluidity emerges 
continuously. For small r we see that Um — —r/u is es- 
sentially proportional to the molecular chemical potential 
/i„j = — r/2, vanishing along r = 0. 

With increasing H, as the TP is approached, u (pro- 
portional to the molecular scattering length) approaches 
zero, so that the onset of superfluidity at the continu- 
ous transition becomes progressively steeper. The point 
r = u = defines the TP (Tt, Ht). For T < Tt {u < 0), 
the superfluid to normal transition is flrst order, with n^ 
and the magnetization M jumping discontinously at the 
phase boundary r = 3u^/16w, yielding phase separation 
at flxed T and M below Tt (Fig. [TJd). 

Close to the TP, however, fluctuations of yield im- 
portant modifications to the mean-field picture, best cap- 
tured via an RG analysis, in the form of logarithmic cor- 
rections to scali ng l29l . l30l . l3ll . |33 | (already observed in a 
different setting [33|). For example, consider the molecu- 
lar density along a line intersecting the TP (i.e., H = Ht 
and u = 0). Along such a line, mean- field theory pre- 
dicts Urn — \/ —T jv, which, along with r (x T — Tt, yields 
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with a logarithmic correction, in contrast to the behavior 
of the standard superfluid transition [35|, described by 
a small exponent a ~ —0.01. In the normal state, for 
i > 0, we find C oc i~^/^, with no logarithmic correction. 
We also find that the divergence of the correlation length 
near the phase transition (recently measured in a bosonic 
cold-atom system [3a|), ^ oc \t\~'^ , reflects the correlation 
length exponent j^ = 1/2 with no logarithmic correction. 
Tricritical scaling also implies a universal shape to 
the coexistence curve near Tt, asymptotically equal to 
the jump in the magnetization SM across the first-order 
phase boundary (recently studied near the TP by the 
MIT group as a jump in the magnetization vs. ra- 
dius [13); we find (Fig. [T]d): 



FIG. 3: (Color online) Plots of the molecular density Um vs. 
position X (normalized to reach unity at the edge), within 
the local density approximation [by replacing r — > (1 — x^) in 
Eq. (|33p ]. for v — 1 with u — (tricritical case, solid curve) 
and w = 1 (above the tricritical point, dashed curve). For 
comparison, the dot-dashed curve shows a standard parabolic 
Thomas- Fermi profile Um oc (1 — a;^). 

Um OC \t\'^^ with t = ^ * and the tricritical exponent 
f3 ~ 1/4, in contrast to /? = 1/2 for the standard su- 
perfluid transition. In fact, the predicted onset is even 
steeper when fluctuations are accounted for, with 
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Thus, fluctuations do not alter the mean-field exponent 
but instead provide a logarithmic factor; this occurs be- 
cause the tricritical upper critical dimension equals the 
physical dimension d = 3 [H, [H [U, [H [H . 

In a cold atom experiment, the tricritical exponent 
/? = 1/4 (and the logarithmic correction) should also 
be reflected in the spatial dependence of the pair den- 
sity in a harmonic trap which is captured, within the 
local density approximation (LDA), via the replacement 
A*m —^ Mm — ^mbfl^x^ (with n the trap frequency). In 
a usual superfluid this yields the well-known parabolic 
Thomas- Fermi (TF) profile n„i oc (1 — x^ /x^p) of a 
trapped gas (refiecting the usual exponent j3 — 1/2), 
with xtf the TF radius. In the present case of a polar- 
ized Fermi superfiuid near the tricritical point, we find 
(setting Xtf = 1 for simplicity): 



,(a;)~ ^/l-x^^\n^''^{l^x^)- 



(6) 



a distinctly different cloud shape (including a divergent 
slope at X -^ 1) for polarized superfluid Fermi gases, as 
shown in Fig. [S] The enhanced symmetry near the TP 
is also reflected in the heat capacity of this interacting 
Fermi gas [S^]; by considering the free energy for i < 
along with the logarithmic corrections, we find 



Coc 1^1-1/2 lni/2|i| 



(7) 



for the heat capacity below Tt, again reflecting the mean- 
fleld tricritical heat capacity exponent a = 1/2 along 



The preceding expressions apply for uniform polarized 
Fermi gases close to the TP (and trapped Fermi gases, 
within the LDA, that have part of their system locally at 
the TP.) Away from the TP, however, such observables 
will cross over to critical (for H < Ht) or first order (for 
H > Ht) behavior asymptotically close to the transition. 



III. MICROSCOPIC MODEL 

To derive these results, we begin with the stan- 
dard one-channel model of two-component resonantly- 
interacting fermions Cj,^ (with cr =| , |): 



^ ^ Z^'^fefCkCTCkCT + ^Z^^ktSi'^k+qiS-qT' 



(9) 
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where S^ka- — ^k — fJ-a-, Cfc — k"^ /2m, and m is the fermion 
mass. The population imbalance is induced by the dif- 
ference in chemical potentials /Zo- = fj, ± H, and the at- 
tractive Feshbach resonance interaction is captured with 
the coupling constant A < 0, the magnitude of which 
increases with decreasing Feshbach resonance detuning, 
and which is best characterized by its connection 
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(10) 



to the scattering length a^. 

We now proceed to focus on the TP of Eq. ^ at fi- 
nite temperature and polarization. We do this by first 
mapping Eq. Q onto the general tricritical free energy 
Eq. (m thereby deriving the coeffcients r, u, and v. 
This involves making the BCS mean-field approximation 
A(x) = A(C|(x)c|(x)), expanding to leading order in the 
magnitude and spatial gradients of A(x), and tracing 
over the fermionic degrees of freedom in the partition 
function following standard methods [33| (generalized to 
H y^ 0). The inaccuracy of the BCS mean-fleld approx- 
imation near unitarity means our predictions for the lo- 
cation of the tricritical point, and the precise forms of r. 



u, and V are not quantitatively trustworthy; however, we 
expect them to be quahtatively vahd (and, they must be 
if a tricritical point indeed occurs). Moreover, our prin- 
cipal interest concerns the power laws and logarithmic 
corrections near the tricritical point, that are indepen- 
dent of the BCS approximation. 

With these caveats we begin by first assuming uniform 
A; the gradient term in Eq. ([T|) will be derived below. 
With uniform A, we can easily trace over the fermion 
degrees of freedom to obtain the following mean- field free- 
energy: 



F = 
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(11) 



-rEEln(l + e-^(^'=+^«)), 



where Ek = yj(f. + jAp. Using Eq. HO]), Taylor expand- 
ing order by order in A, and evaluating the momentum 
integrals yields: 



i^ = ^o + ^2|A|2 + iF4|A| 
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(12) 

(13) 
(14) 
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with ipn {x) the polygamma function and Re denoting the 

real part. In evaluating these integrals, we have assumed 

/i > and expanded the fermion density of states to 

leading order at the Fermi surface (yielding the factors 

3/2 1/2 
of Po = ™/K tt.3 ) 1 ^-n approximation that breaks down in 

the deep BEC regime, where ^ becomes negative. 

The tricritical point occurs when V2 = V4 = 0. Ex- 
amining Eq. (|14p . we see that mean- field theory pre- 
dicts a universal ratio between the chemical difference 
Hi and temperature Tt at the tricritical point, given by 
Re[^4(i + ^)]=0'Or, 



Tt 



(16) 



valid at any detuning (within the preceding mean-field 
assumptions). Equation ()16|) can be combined with 
Eq. (fT3|) to obtain a prediction for the ratio Tt/ ^jl as a 
function of the scattering length as that gives 






0.35, 



(17) 



at the unitarity point a^^ =0. 

Our next task is to obtain the gradient term in the 
free energy. We do this directly following the textbook 



derivation of the GL free energy (as discussed in Ref. ^37), 
that essentially allows spatial variations in the quadratic 
coefficient. Here, the only difference is the nonzero chem- 
ical potential difference H . We obtain an iJ-dependent 
coeffcient of the gradient term that, taking H to be given 
by its value at the tricritical point Eq. (fT6|) . yields 



F = d^R 



7m^i;^C(3)c^ 
967r4T2 



ivAr 
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+ 2^4|A|4 



V2\/^\ 



3^6|A|' 



(18) 



where Cm — 0.0789 and we defined the Fermi velocity 
v-p = 27r^po/m^. 

To reduce F to the form of Eq. ^, we define A(x) = 

riip{x.) with r] ~ 3.3( ^ Jrs) ) chosen to fix the coefR- 
cient of the gradient term. Then we'll have r = 2V2rf' for 
the quadratic coefhcient, u = 2V4rf for the quartic coeffi- 
cient, and V ~ 2VQrf for the sixth order coefficient. Near 
the tricritical point at unitarity, these are approximately 
given by: 



2po?7 (^ + 0.86/i), 
-0.16^/z, 



.0064 
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(19a) 
(19b) 

(19c) 



and Po = 



is the density of 



where h = ^^ „... ^^ - -^j^^^ 
states at /i. Thus, as already mentioned, r vanishes along 
a line in the T-H plane, (that is the continuous super- 
fluid to nonsuperfluid transition occuring for T > Tt) 
and u vanishes at H ^ Ht- The sixth-order coefficient is 
approximately constant near the TP. 



IV. RENORMALIZATION GROUP 

Having provided an approximate connection between 
the microscopic one-channel model Eq. ^ and the free- 
energy functional Eq. ([1]), we now proceed to analyze the 
enhanced critical fluctuations near the tricritical point of 
F. Such fluctuations, and the concomitant logarithmic 
corrections to scaling [2^, [30, [3l|j [sj , are best analyzed 
using the renormalization group (RG), which incorpo- 
rates the effect of fluctuations neglected in the mean-field 
approximation. It is convenient for such an analysis to 
set h /mi, = 1 (equivalent to measuring lengths in dif- 
ferent units), so that the gradient term of Eq. ^ has 
coefficient 1/2: 



F = I d^x[i| VV'P + \r\i^? + \u\^^ + \v\i^\'] . (20) 

We first note that the model Eq. ([^0)1 is defined for mo- 
menta below an upper cutoff A, i.e., it is coarse-grained 



on length scales larger than ^ %/ K. (For a unitary po- 
larized Fermi gas, we expect A ^ /cpi with fcp the Fermi 
wavevector, on dimensional grounds.) The perturbative 
RG proceeds by integrating out states close to A. Thus, 
we split the Fourier transform V'(p) — J d^xe^'^'^ip{x.) 
into low and high momentum modes ipip) = '0<(p) + 
■0>(p) with ^/'<(p) defined for < p < A/6 and ^/'>(p) 
defined for the shell of momenta A/b < p < A with 6 > 1. 
By evaluating the trace over ip^ (x) , perturbatively in u 
and V (focusing on the crucial In A-divergent terms), we 
derive an effective theory for ip^ {x) that is of the form of 
the original model but with renormalized couplings. 

Let's illustrate this procedure for the contribution to 
the effective Hamiltonian from expanding the sixth-order 
term to 0{v^). This is: 
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d^Xid''x2{\iJ<{x,)+iP^{x,)\^ (21) 
x|V'<(a;2)+^>(a;2)|^)>, 



where the subscript > indicates the trace over the high- 
momentum modes. This trace is evaluated with the help 
of the two-point Green function 

{^l{xi)^>ix2)) - 2G>(a:i - X2), (22) 



which has the Fourier transform 



G>(p) = 



1 



p2 _|_ J, 



(23) 



for A/b < p < A. This yields: 



F^2 = -224— / d^x\i;<{x)f / d^x'G>{x'f, (24) 

so that the O(w^) term clearly renormalizes the bare 
sixth-order coupling v. To proceed, we must evaluate 
the final integral in Eq. (j24p . If the momenta appearing 
in G{p) were unrestricted, then this integral would be 
In A divergent, and given by: 
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(27r)3 7 {2TT)^pl + rpl + r{pi+p2)^+r' 

(25) 
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where we took r > for simplicity. In the present case, all 
momenta are restricted to the high-momentum window 
A/6 < p < A; for this case wc find 



d-'x'Gyix'y c^I{A)-I{A/b) 
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IGtt' 



■In 6, (26) 



thus yielding 

F„ 
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ln6— / d^x\i;<^{x)\'^. 



(27) 



To complete the RG procedure, we must define new mo- 
menta p' = bp to restore the original cutoff A, equivalent 



to defining x' = b~^x. As usual, this also causes a rescal- 
ing of the kinetic energy term in Eq. (|20p that must be 
absorbed into a new fermion field tpi^')- 



ij<{bx') =6-3^(x'). 



(28) 



With these manipulations, the final effective Hamiltonian 
is exactly of the form of Eq. (pO|l but with the renormal- 
ized coupling 



v' = V — cv In 6, 



(29) 



with the numerical coefficient c — ■^. Upon iterating the 
RG procedure, and including similar renormalizations for 
the quadratic and quartic coefficients, we find the RG 
equations: 

dv{b) 
dlnb 
du{b) 

dlnb 
drjb) 
dlnb 

for the running coupling constants u(6), v{b) and r(6), 

with C4 = -% and C2 — -J-^. These RG equations, consis- 

I — I if^ 
tent with previous results reported in Refs. |31|, |3J, can 

be integrated to yield 



= -cMbf. 


(30a) 


= U{b)[l - CiV{b)], 


(30b) 


= 2r{b) - C2u{bf , 


(30c) 
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v{\ + culn6) , 
;ln6) 
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(31a) 
(31b) 



lli2 



r(6) = 62(r-t--— [l-(l + CTln6)^]), (31c) 

with the latter equation being given by rib) ~ 6^r close 
to the tricritical point. 



V. MOLECULAR DENSITY 

Having computed the RG equations, our next task 
is to combine these with scaling relations for various 
experimentally-observable quantities. We begin with the 
molecular density [2^ n^ = IV'P: which satisfies the scal- 
ing equation 

nm(r, u, v) = b~'^nm{r{b),u{b),v{b)), (32) 



following from Eq. (I28p . Below, we use the shorthand 
nm = b~^nmR for such an equation. The left side of 
Eq. ((32)) is the physical molecular density, while the right 
side is the density in the renormalized system. The RG 
strategy is quite simple: Although mean-field theory is 
invalid close to criticality, we can choose 6 = 6* such 
that the renormalized system is far from criticality where 
mean-field theory is accurate. This is seen directly from 
Eq. (I3ip : with increasing 6, |m(6)| and |r(6)| grow large 
while i;(6) — > (validating perturbation theory). In the 



critical and tricritical regimes, it is sufficent to take the 
condition |r(6*)| ~ 1 which yields 6* ~ l/|r|^/^ for the 
RG condition. Then, using the mean-field result [i.e., 
Eq. (g])] for the right side of Eq. ([321), we find for the 
molecular density: 



~2^ 



where 



A\r\v 



l + culn|r| 2, 



(33) 



(34) 



is a logarithmic correction to scaling. Equation ([33]) de- 
scribes the molecular density near the tricritical point 
along lines that intersect the critical line or the tricritical 
point. In the tricritical regime H ^ Ht, where u ^ 0, 
this yields Eq. ([5|) for the leading temperature depen- 
dence. As noted in Sec. HH Eq. ([33]) imphes a distinctly 
different spatial dependence of the molecular density in 
a harmonic trap, seen in Fig. [31 with a spatial variation 
described by Eq. ([6]). 

For any T > T^, above from the tricritical point, 
Eq. ([33]) crosses over asymptotically to critical behavior 



l«3/7 



(35) 



therefore yielding a standard TF profile, with Tim ex; 
(1 — x'^ /x'^p) near the edge. However, close to the TF, 
this crossover will only occur extremely close to the edge 
of the cloud. As the tricritical point is approached from 
the critical regime m > 0, the slope of the boundary is 
predicted to diverge as m -^ according to Eq. (I35|) . as 
seen in FigJ31 This vanishing of u can be interpreted as a 
vanishing [l6j of the effective molecule-molecule scatter- 
ing length a^n (recently measured [3a|); using Eq. (j31b[) 
along with the above renormalization condition gives, at 
T = Tt and for h ^ 0, 



a„,oc\h\ln'^^^{\h\~') 



(36) 



for the vanishing of the molecular scattering length close 
to the TP. 

We now turn to the behavior of the molecular den- 
sity across the phase transition in the first-order regime 
(T < Tt) near the TP. As noted above, in the continuous 
regime near the TP, the superfluid order parameter ex- 
hibits a sharp onset at the phase transition, resulting in 
the steep edge of the cloud shape shown in Fig [3] In the 
first-order regime, the order parameter jumps discontin- 
uously across the transition, resulting in a jump in the 
molecular density at the edge of the molecular cloud. 

The molecular density in the first-order regime can 
similarly be computed using the RG; in this case, how- 
ever, we must choose a more generally- valid renormaliza- 
tion condition ^39)]. We'll take this to be when the cur- 
vature of the free-energy minimum -yh^ in the ordered 
state is large (i.e., order unity); this yields 



1 



1 



4r»v* — u^\/ul ~ 4r»f*] , (37) 



for the most general condition in the ordered state, with 
v^, = ^(6*), etc. Close to the continuous transition or to 
the tricritical point, Eq. ([37]) for &* is essentially equiv- 
alent to 6* ~ l/jrl'^/^ as used above. In the first-order 
regime m < 0, however, Eq. ([37]) is, approximately, 

1^2ul/v,, (38) 

giving our renormalization condition in the first-order 
regime. To compute the jump in the order parameter 
or molecular density across the phase boundary, we need 
use this condition to find the first-order phase boundary, 

o 2 

which in mean- field theory is given hy r — -^^. Com- 
bining this with Eq. ([35)1 then yields 



(39) 



r = ^£* 

for the renormalized first-order phase boundary. Along 
with Eq. (|32[) . we can then obtain the jump 



Srir. 



3H 
4u 



1 + cv\n 



V2h 



(40) 



in the molecular density across the first-order phase 
boundary, which vanishes as u ^ at the TP as ex- 
pected. 



VI. MAGNETIZATION 

We now turn to the local polarization, or magneti- 
zation, M = n| — Tij^, a crucial observable in present- 
day experiments, that is related to the free energy via 

M = 



OF 
dH- 



In deriving Eq. ([!]), we dropped the overall 
constant term Fq in the free energy; we must now rein- 
state this term which yields a fermion-only contribution 
to the magnetization that we denote Alf. Thus, we have 



M = M 



1 dr 



1 du 



1 dv 



(41) 



^ 2ai/'^' 4ai/'^' eaij'^' ' 

for the magnetization. 

Near the tricritical point where ji/'l vanishes, Mf will 
be the largest of the terms in Eq. ([4T|) but not show any 
significant signature of the phase transition. However, 
one observable that will show a sharp signature across 
the phase transition is the jump in the magnetization 
across the phase boundary in the first-order regime. This 
quantity, that is directly measurable as the jump in mag- 
netization as a function of radius in a trapped polarized 
Fermi gas [Tjl, is equivalent to the width of the coexis- 
tence region SM = Mc2 — Md below the tricritical point, 
see Fig. [Hd. 

The O(l^p) term provides the dominant contribution 
to 6M. Using Eq. ([40]) for the jump in \i/j\'^ across the 
phase boundary along with Eq. ([35]) . we obtain 



SM 



a/3 dr 
TdH 



1 



) Ini 



1/2)1/2^ 



(42) 



for the jump in the magnetization in the first-order 
regime near the TP. 



VII. HEAT CAPACITY 

In the present section, we use the RG to obtain correc- 
tions to mean-field theory for the free energy and heat 
capacity C = —T^^. Cfose to the phase transition, we 

-1)2 771 

can take C on ^^ , since the leading T dependence is via 
r. Using the free-energy scaling relation F = 6^^Fr, we 
find 

C{r,u,v) ^bC{r{b),u{h),v{h)), 



(43) 



for the RG equation for the heat capacity. We begin with 
the heat capacity in the ordered phase. In this regime, 
we'll use the mean-field result for the heat capacity for 
the right side of Eq. (|43|) . which is obtained by differen- 
tiating the mean-field free energy Eq. ^ and using the 
stationarity condition Eq. ([3]): 









(44) 



with the angle brackets reflecting the thermodynamic av- 
erage with respect to F. We can obtain the leading-order 
approximation to (IV'P) using the normal-state Green 
function: 



d^k 



f27r)3 P 



(49) 



which, upon differentiating, gives for the heat capacity 
in the normal state: 



c-i-.(l;f ' 



'dT^ SttV?' 



(50) 



Using this result for the right side of Eq. (|43p . and taking 
r(6*) to be given by Eq. (|31cp . we obtain for the heat 
capacity in the normal state (r > 0): 



C~T, 



dr s2 1 



dT' 8tt 



1/7n 



(51) 



+ i^a-er") 



Taking the RG condition 6* = l/-\/|r| (appropriate for 
the critical and tricritical regimes, on which we shall fo- 
cus), we finally obtain in the ordered state r < 0: 



1 (jr 9 



1 



(45) 






Av\r\er 



for the fluctuation contribution to the heat capacity 
which, for H = Ht on a, line crossing the TP, gives 



which, along a line intersecting the tricritical point (u = 
0), immediately reverts to Eq. (j50p . i.e., we find no loga- 
rithmic corrections in the normal state at the tricritical 
point so that, for H = Ht, 



C< 



1 



(52) 



in the normal (nonsuperfluid) phase. 



^^^'~^df^ 



dr -.2 yirijrp^ 



(46) 



VIII. CORRELATION LENGTH 



that reduces to the previously quoted formula for the 
tricritical heat capacity, Eq. ([7|, once we use r ex i [as 
follows from Eq. (|19ap ] 

In the critical region above the TP (for T > Tt), 

Eq. (gni) gives [13 C oc (2u)-i£?^^. In the absence of the 
log correction, this would represent the usual mean-field 
specific heat jump at a superfluid transition; however, 
the true asymptotic behavior in this regime is known to 
reflect a very small heat capacity exponent a |35j . 

In the disordered (normal) state, the mean-field free 
energy vanishes (since ip = there). Thus, to obtain a 
nonzero result for C at r > we must go beyond mean- 
field theory. The free-energy in the normal state, !F, is 
given by a functional integral over the field tp 



T = -Tin / D^e 



-FjT 



(47) 



with the action F given in Eq. (pOl) . Near the transi- 
tion, we can assume that the leading T dependence comes 
through r, which yields 



C 



■-€n>'^ 



(48) 



Finally, we consider the divergence of the correlation 
length ^ near the phase transition, governed by the decay 
of the superfiuid correlation function at large distances. 
Here we shall focus only on a line crossing the tricritical 
point, i.e., we take u = 0. In the normal phase, we have 



(V'^(a;)V(0))=2G(x) 



(53) 



with the Green function satisfying the scaling relation 
G{x) = b^^GRib-^x) [as follows from Eq. ^]. Using 
the perturbative result for Gr{x) and the renormaliza- 
tion condition b^ = 1/ \/r gives 



(Vn^WO))=^e-v^l-l, 



(54) 



identical to the perturbative result [due to the trivial 
scaling of r near the TP Eq. picp ]. which, upon examin- 
ing the argument of the exponential gives the correlation 
length ^ oc r~^''^, and the exponent v = ^. A similar 
analysis in the ordered phase (expanding the free-energy 
around the mean-field solution) gives ^ ex |r|~^/^ below 
the TP. 
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IX. CONCLUDING REMARKS 

To conclude, we have computed various experimen- 
tal predictions for the behavior of a polarized superfluid 
Fermi gas near its tricritical point using a GL functional 
that is generally valid near the tricritical point even in the 
unitary regime. We presented results for numerous ob- 
servables in cold-atom experiments, including the onset 
of the superfluid order parameter at the transition, the 
magnetization jump across the first-order phase bound- 
ary, the heat capacity, and the correlation length. In 
contrast to the standard superfluid transition, in which 
critical exponents deviate from their mean-field values 



(exhibiting anomalous values), for tricritical points the 
mean-field exponents are predicted to be exact, but with 
logarithmic corrections. These predictions should pro- 
vide sharp signatures of the tricritical point in polarized 
Fermi gases. 
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